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The Two- Dimensional Motion of a Plane Lamina in a Resisting 

Medium. 

By S. Brodetsky, M.A., Ph.D. 
(Communicated by Prof. G. H. Bryan, F.E.S. Received July 20, 1916.) 

The General Prohlem. 

1. A plane lamina, AB (fig. 1), is perpendicular to the plane of the paper, 
and moves in two dimensions in this plane. We use the following 
notation : — 

V is the velocity of the centre of mass, C, of the lamina ; 

u, V are the velocity components along and perpendicular to the lamina; 

6 is the inclination of the lamina to some fixed direction, CD ; 

a is the angle between the velocity V and the direction of the lamina — 
the '' angle of attack '^ ; 




Fig. 1. 



(j) is the angle between the velocity Y and the fixed direction, CD, so that 

E is the normal air resistance, supposed to act at C ; 

G- is the couple which, with R, represents the effect of the air resistance, 
tending to diminish ; we measure forces in dynamical units, i.e., in poundals 
or dynes ; 

M is the mass of the lamina and its appendages, so that — 

Mg is the weight; 

I is the moment of inertia of the lamina, etc., about an axis through C 
perpendicular to the plane of motion ; 
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X, Y are the components of the external applied forces, as gravity or 
propeller thrust, along the direction CD^ and perpendicular to CD down- 
waids, acting at the centre of mass. 

We suppose the centre of mass to be in the lamina at C, and the centre of 
pressure for a = 90*^ to coincide with the centre of mass when there is no 
rotation, so that in this standard case the couple G is zero. 

The equations of motion can be written in either of the following forms :— 
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It is sometimes more convenient to use the form (I), and at other times 
the form (2). 

2. The quantities R and G are complicated functions of the size and 
shape and disposition of the lamina, of the angle of attack, a, of the 
velocity V, and of the rotation dOjdt. Experiment shows that air resistances 
vary as the squares of the rates of displacement. It follows that we ought 
to use the general formulae 



E = vy.(4f) 



G = vy.(a,l|), 



(3) 



where /i, /2 are functions whose coefficients involve the geometrical circum- 
stances of the lamina. Such general forms would, however, render the 
equations of motion quite unintegrable. We shall |;herefore use only the 
first one or two terms in the functions supposed expanded in Fourier's series. 
If the lamina is symmetrical about a line through C perpendicular to the 
plane of motion, we may write 



E = KSV2SA„sin(27i+l)a + 2KSV^SA„'cos(2yH-l)a 

^^0 



00 



M^„ . . „ /dd^^ 
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G = KSV2SE„sin2«a-|-2KSV^SE„'sin2?i«+KSO(-,, 
1 di I \dt 
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where S is the area of the lamina, K is a resistance constant depending upon 
the disposition of the area and the nature of the medium, and the A's, E's, 
and C are numbers. Experiment shows that in each infinite series only the 
first few terms are at all important. 

But we can simplify our problem if we suppose, as is generally the case 
in actual aeroplanes, that the dimension, AB, of the lamina is small, com- 
pared to the dimensions of the appendages of the lamina. The effect of this 
supposition is that we may altogether neglect the terms in the air 
resistance that depend upon powers of the dimension, AB, higher than 
the first or second. If then, the angular velocity is not large, we may put 
G zero, and make E proportional to V^. In this case of a narrow lamina 
we shall in the first instance use the sine law 

KSV^Ao sin u (6) 

as a hypothetical value for E. But, in order to make our results more 
applicable to practice, we shall consider the composite form 

E = KSV^Ao sin «, from « = 0° to « = 30^, -i 

\ (6) 

E = JKSV^Ao, from « = 30^ to a = 90^, J 

with the obvious changes in the sign of E when a is in other quadrants. 
This agrees tolerably closely with the measurements of Eiffel and others, as 
well as with the formula suggested by Duchemin. 

Finally, in accordance with a suggestion by Prof. Bryan, we shall investigate 
the problem with 

E = f KSV^Ao sin a cos «, from « = 0^ to « = 65^ 42', 

E = JKSV^Ao, from « = 65° 42' to 

This also agrees tolerably with experiment in the case of a square lamina, 
and has the advantage of simplifying the analysis (see § 8). 

If the lamina is not quite so narrow, we have to take account of the shift 
of the centre of pressure. If we still suppose the rotation to be moderate, 
we can use for E any of the values (5), (6), or (7). But for G- we write 

G = KSV2Eisin2a, (8) 

which will be found to agree fairly well with the measurements of Soreau 
and others, for all values of the angle of attack a. 

Lamina with Large Moment of Inertia, Motion under No Forces, 

3. We first consider the motion of the lamina under the sole influence of 
the air resistance, there being no other external applied forces. 

Very Narroio Lamina. — If the lamina is very narrow, whilst the moment 



65^42', ^ 

I (7) 
a = 90°. J 
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of inertia, I, is considerable, we are justified in assuming that the resistance 
couple has no appreciable effect. Using the sine law for E, the equations of 
motion (2) give 

dt dt 

dv , dO KSAot72 • I /f^\ 

^'^ = 0. 



dt 

dO 
We thus get — - = w, a constant, 

dt 



J 



so that the first two equations of (9) become 

d% 

j> (10) 



dt 



---+a}U + Av(ti^-i-V'^)i — 0,J 

where A is the constant KSAo/M. 

Equations (10) give ~~(u^ + v^)^ == v, (11) 

d% CO 

a relation between it and v. This can be written in the form 

r/V 

since V^ = t{?'\-v^, 

h is the constant A/w. 

The differential equation (12) can be used to find out the nature of the 
motion. In spite of its very simple form, it does not appear to be integrable 
in exact terms. We can, however, use it to obtain a graphical relation 
between n and v. 

It is convenient to adopt the convention indicated in fig. 2. We take « 
to be positive always, Further, u is taken to be positive or negative 
according as it is from C to B or from to A. Finally, we call v positive |f 
it is on the side B to A of the rotation, and negative if it is on the side A 
to B. The effect of this convention is that k can be taken to be positive in 
all cases. 

This constant, h, in the equation (12) only affects the dimensions of the 
curves represented by this equation, and not the shapes. Thus, if we plot 
the curves represented by the equation 

y-=~^^ (13) 

du 
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Fig. 2. 

we have only to multiply each radius vector by the constant Ifh, and we 
obtain the curves (12). The (%v) curves (13) constitute a singly infinite 
family, no two members of which can pass through the same point. The 
family as a whole possesses symmetry about the origin, for the equation (11) 
is unaltered if we change the signs of both u and v, 

4. In order to plot the curves (13) it is useful to find the locus of the 
points of inflexion on the curves, as well as the locus of the points at which 
the curves are parallel to the axis of u. From (13) we get after some reduction 

Hence the points of inflection lie on the curve 

^4_^,^^V-Y^ = 0. (14) 

This locus is plotted in fig. 3, dots and dashes. 

The points at which the (% v) curves are parallel to the u axis are given by 
the condition 

dv . rx 

d% 

■i.e., tb-^vY = 0, (15) 

bv (11). This locus is also traced in fig. 3, broken line. It is asymptotic to 
the locus of inflexions, both loci being asymptotic to the lines v =: ±1 as 
shown in the figure. 

Some of the {%,v) curves (13) are shown in fig. 3. Each one consists of a 
portion beyond the locus of inflexions, convex to the u axis; then there is a 
point of inflexion, and after that each curve assumes the form of a spiral, 
everywhere concave to the origin, and reaching the origin after an infinite 
number of convolutions, whose ultimate shape is circular. The limiting 
case, in which the point of inflexion is at infinity, as well as the furthest 
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intersection with the % axis, is represented by the two dotted curves in the 
figure. These are asymptotic to the locus of inflexions and to the locus of 
maxima and minima — points at which the curves are parallel to the tt axis — 
and are situated between these loci. 

To adapt the figure to the more general value k in (12), we let the distances 
of the asymptotes from the % axis represent the number IJh, i.e., w/A, the 
scale of u being altered in the same ratio. 

5, It is readily seen that the progress of changes in u and v obtained by 




Fig. 3. 



moving along any {u,v) curve inwards to the origin, corresponds to the 
continuous progress of time forwards. Thus we have by (10) 



dt 



= (aVj 



where according to our convention m is positive. Hence when v is positive, 
u increases with t, and when v is negative, u decreases as the time progresses 
forwards. Eeference to fig. 3, shows that these facts imply a continuous 
approach to the origin along any curve as t increases. 

If we take any point {%t,v) there is a definite {u,v) curve passing through 
it. Thus, with given initial values of the velocity components we get to know 
exactly the changes that take place in the components, by following up this 
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curve till we come to the origin. The radius vector to any point in the (u, v) 
figure represents the resultant velocity, V, and the angle that the radius 
vector makes with the u axis is the angle of attack a at any instant. We 
see then that the velocity, Y, continuously decreases to zero, whilst the angle 
of attack passes through cycles of progressive changes 27r to 0. 

It only remains to introduce the time element, and thus also incidentally 
the angle 0. By means of the relation 

dt = ~ , (16) 

0)V 



given by (10), we can follow out the changes in the various quantities u, v, Y, a, 
corresponding to gradual progression of the time, and therefore also of 0^ 
since d0/dt = co, 

6. We have, therefore, developed a method which can be applied to plot 
the path of the centre of the lamina, as well as the direction of the lamina, 
with tolerable ease. We reserve the carrying out of the process to a later 
stage, where we shall also take into account the effect of gravity. But we 
can obtain some information about the path of the centre of mass by means 
of the following reasoning. 

In the second equation of (1), put 

X = Y = 0, R^AMY^sina. 

We get -^ = —AY sin a cos 5i, 

^ dt 

so that -^ = — A sin a cos a, 

ds 

where ds is an element of length of the path. Plence 

Asm 2a ' "^ 

where p is the radius of curvature of the path. 

Again, since 

d(h d(h dt 1 dd) A 

j_ ^^ X rzz ' —^ ,— cos OL 

du dt dtt (ov dt CO 

the increase of (j) between any two values Ui, U2, is given by 

A r''3 

(j)2--<j>i ~ QO^adtk (18) 

From (17) it follows that the radius of curvature of the path is never less 
than 2/A numerically, it being quite unaffected by the progressive decrease 
in Y. Also, whenever a is a multiple of a right angle, and this happens an 
indefinite number of times before the motion of translation is destroyed, there 
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is a point of inflexion, p = oo . Again, from (18), we deduce that the angle 
<^ oscillates with diminishing amplitude about some limiting value. We thus 
arrive at the following general conclusions : — 

If the initial velocity is only moderate, we get an oscillating path, and the 
oscillations gradually die down. If the initial velocity is big enough, we get 
one or more series of loops, and then the diminishing oscillation. 

Meanwhile, in all cases, the velocity diminishes continuously, and the 
angular velocity remains constant. Hence, ultimately, the rotary part of the 
motion predominates, until finally, after an infinite time, only the constant 
rotation is left. This is, of course, what would be expected from the principle 
of energy. 

7. The sine law for the normal air resistance is inapplicable in practice for 
values of | sin a | greater than J. For larger values of | sin a | , i.e., for the 
parts of the (2^, 'y)-plane where \v/u\ > 1/^/(3), we need to adopt some other 
hypothesis, a suitable one being 

E = JKSV^Ao = JAMV^ (19) 

using the notation of § 3. When v is positive, (2) become 



du '^ 

dt == "^'^ 



^^ + mt + ^A(u^ + v^)= 0, 



y (20) 



dt 

ft) being the constant angular velocity. The solution for v negative is 
obtained by symmetry about the origin. 
We get for positive v the equation 

dtt CO 

i.e., V2 = a-A^/-, 

where C is an arbitrary constant. This can be written 

kY = Ge~'-^^ (21) 

k being the quantity A/co as in § 3. We may, once more, put k = 1 in the 
first instance, and then adapt the resulting curves to the general value of k by 
multiplying their dimensions by co/A. Thus, in the part of the upper half of 
the (7^, 'y)-plane between the lines v/to = + 1/\/(3), we plot the curves 

Y = C^-i^. (22) 

When V is negative, i.e., in the lower half of the {u, v)-plane between the 

lines v/u= ±1/^(3), we use 

V = Cei\ (23) 

In the remainder of the plane we use the curves given by § 4, plotted in 
VOL. xc V. — A. 2 T 
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fig. 3. We thus obtain a ^et of composite {u, v) curves, which represent 
with reasonable accuracy the motion of the narrow lamina under no forces 
but subject to the resistance of the air, or of any other medium. As in | 5, 
we have ih, v, V, a in terms of the time, whilst the direction of the lamina 
is given by dOfdt = co. 

The composite curves are plotted in fig, 4. The points of inflexion lie on 
the curves marked in dots and dashes, and the points wdiere the curves 
are parallel to the u axis lie on the curves marked in dashes. The dotted 
curves represent the limiting cases, when ie starts off at + co. 




Fig. 4. 



8. For a square plate, or one in which the aspect ratio is a little greater 
than unity, Prof. Bryan suggests the use of the law of resistance, E, varying 
as the sine of twice the angle of attack, for values of this angle up to some- 
thing greater than 45^, in accordance with the measurements by Dines. 
Although the maximum in Dines' curve (see Lanchester's ''Aerodynamics/* 
fig. 101) is about 35° still the hypothesis 
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E = 4 AMV^ sin ol cos a, from <% = O'' to a = 65° 42', -^ 

I (7) 

= 1 AMY2 . from a = 65° 42' to ^ = 90°, J 

with corresponding values in other quadrants, is a fair representation, and 
can be used to examine the effect of the maximum. The maximum is 
about |- times the more or less constant value alter the drop. Also 
I sin 65° 42' cos 65° 42: is |. 

With the sine twice angle hypothesis, equations (9) give 



du -^ 



> (24) 

dt ^ -^ 

When It is positive, equations (24) give 

:7-+-~+t;"^^ = 0, (25) 

du V 6 0) 

which on integration gives 

^lm-'\og{l + ^hv)-\-l{UMf = const, (26) 

k being once more the constant A/cy. We again put A; = 1, and then adapt 
the results to all values of k by multiplying the dimensions of the {u, v) 
curves thus obtained by the quantity Ijk. 

Hence, in the right half of the {u, 'y)-plane, between the line? 
vju = Hhtan 65° 42', we plot the curves, 

|^-log(l + |'y) + i(I^O^ = const. (27) 

For u negative, i.e., in the left half of the {u, 'y)-plane, between the lines 
q^ju z=z +tan 65° 42', we use symmetry about the origin, or plot 

-|'y-log(l-|'y)+K-i-^)^ = const. (28) 

It is important to notice that the constants in the equations (27-8) 
may also have complex values, so that we get (%i, v) curves covering the 
whole of the space to which the sine twice angle law applies. 

In the remainder of the (u, i^j-plane, i.e., above and below the lines 
v/u = ± tan 65° 42', we use the curves given by § 7. We thus get once 
more a set of composite (u, v) curves, which can be used to discuss the 
motion. These are plotted in fig. 5. The lines in dots and dashes are the 
inflexion loci, and the dashed arcs are the points where the curves are parallel 
to the u axis. The dotted curves are the limiting cases, starting off at 

2 T 2 
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u = -foo,'y= +1 ft)/A. In these limiting cases we have v constant, whilst 
^6 increases from — oo or decreases from + oo , until the value 



— 3 a) 



q:^^ cot 65" 42' 

4 A 



is reached, when v also begins to vary. 




Fig. 5. 

This curious result is obvious from the equations (24). It means that under 
the sine twice angle law, if % = + f co/A, v remains constant. If this law 
were to hold for all values of a, this constant value would be maintained 
till u reached the value zero — not later, as is clear from considerations of 
energy which necessitate a continuous decrease in the resultant velocity. In 
the present investigation, the constancy of v stops at the limit of the applica- 
bility of the sine twice angle law. 
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Wider Lamina. 

9. When the moving lamina is not so narrow that we may altogether 
neglect the shift of the centre of pressure, we must take account of the 
resistance couple. If we still suppose the rotation to be small, we 
can, using the sine law or a substitute, again reduce the equations of 
motion to ordinary differential equations of the first order, similar to 
but slightly more complicated than those already given. These can be 
solved graphically, and several interesting properties of this type of motion 
deduced. 

Motion under any Forces, 

10. It is evident that, if we introduce additional applied forces, the 
equations of motion will be even less tractable than those for motion under 
no forces. Before adapting the graphical method of this paper to the 
more difficult problem, we must mention a method that has often been 
used in the study of particular cases. The method consists in assuming 
some arbitrary but convenient relation between two or more of the variables 
that define the motion. 

Although such an assumption must be quite artificial, and in general not 
quite justified, nevertheless the consideration of such artificially constructed 
problems can lead, and has led, to important results of practical interest and 
application. In connection with the subject of the present investigation, 
reference should be made to a special " artificial " problem worked out by 
Kene de Saussure (Appendix to Langley's "The Internal Work of the 
Wind," ' Smithsonian Contributions,' 1893). De Saussure assumes that the 
angle of attack is inversely proportional to the square of the velocity. This 
amounts to presupposing some agent introducing the necessary couple at 
every instant, so as to give the required relation between a and V. In this 
way de Saussure is able to obtain a mathematical illustration of Langley's 
theory of the soaring of birds. 

The most important case of such an artificial problem is afforded by 
Lanchester's " Phugoid Theory " (' Aerodonetics,' Ch. II., etc.). In order to 
obtain an integrable case of the motion of a plane lamina in air under 
gravity, Lanchester presupposes .applied forces which have the effect of 
supplying the energy that is lost owing to the dissipative effect of the air 
resistance. It is clear that this cannot be realised in practice, yet the 
paths thus obtained — the phugoids — have the great advantage of being near 
approximations to fact when not far removed from the horizontal, and 
hence Lanchester was able to discover the correct condition for the longi- 
tudinal stability of an aeroplane. 



528 Dr. S. Brodetsky. The Two-dimensional 

11. In order to apply the graphical method to the case where there are 
applied forces such as gravity or a propeller thrust, we shall need to make 
use of the curves connecting the velocity components % v. "We shall use 
the composite curves of figs. 4 and 5. These give us what we may call the 
'' fundamental {it, v) curves for a narrow lamina " of large aspect ratio, and 
of aspect ratio about unity, respectively. The fundamental curves refer to 
no forces. It is still assumed that the angular velocity is moderate, so that 
we need not introduce any terms into the resistance force and couple 
involving the rotation. The external applied forces are supposed to pass 
through the centre of mass. 

When there are impressed forces besides the effect of the medium, let us 
consider what happens during an interval of time, dt. During this time the 
representative point in the (w, v) plane travels from its initial position to 
some neighbouring position on the {u, v) curve, passing through the initial 
position. At the same time, lo and v are increased by amounts 

(-COS ^ + X sin 0]dt, / Jcos6>-?:Sin 6\dt, (29) 

\M M / \M M / 

X, Y being the impressed force components, as in fig. 1. We write these 
increments in the form 

Xcos^ + Ysin0 7 Ycos^—Xsin^, .qa\ 

— d% — ^^^^ cm, (oU) 

If, then, in fig. 6 the point Pi represents the motion at the beginning of 



11 



./ 




Fig. 6. 



tlie interval dt, and the curve shown in the figure is the (u, v) curve through 
Pi, during the interval we go to P2, where 

'NT2=^'^dt, 

dt 
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Now measure off P2P3 parallel to the u axis, and P3P4 parallel to the v axis, 
so that 



VV _Xcos^ + Ysin^ ^p ^ 
^ ^ M{du/dt) ' '^ 

pp _Ycos^— Xsin^ -^p 



(31) 



c^^ 




(i'^/' 


rf^ 


= ft), 


f^^ 




P2P4 . 


^f/ 



both quantities being taken with their proper signs, and du/dt and 
referring to the fundamental curve through Pi. The representative point 
in the (u, 'y)~plane will therefore in the interval dt go from Pi to P4. The 
point P4 is on a new fundamenfcal curve. We repeat the process indefinitely. 
In this way we obtain the aditcd (it, v) curves, when X, Y are given at all 
points. 

The important case is where X, Y are due to gravity, i.e., 

X =0, Y = M^, 

so that P2P3 = ^^^i^i^NPs, P3P4 = ^^^ NP2. (32) 

dtf^/dt du/dt 

12. We consider here the narrow lamina of large aspect ratio, using the 
fundamental curves given by fig. 4. In this case we have proved 

Hence -L^^Jl^^ (33) 

and the angle that P2P4 makes with the axis of v^ measured clockwise, is 

= 0Q + a)t, (34) 

$0 being the value of ^ at ^ = 0. 

Using the value of P2P4 in (33) in the direction given by (34), we get the 
actual (ii, v) curves, and at the same time have marked on them points 
corresponding to equal increments of time, and also of equal amounts of 
rotation. It is now a comparatively simple process to plot the actual path 
of the centre of mass of the lamina, as well as the direction of the lamina at 
each instant. 

13. Terminal motion would be defined as one in which the circumstances 
of the motion, namely V, oe, 0, recur periodically. Such a case cannot occur 
if there are no applied forces, for, owing to the resistance of the medium, 
energy is being dissipated continually throughout the whole of the motion. 
But with external applied forces periodic paths can be obtained. A (u, v) 
curve plotted from the fundamental curves in the manner described in 
§§ 11, 21, which, starting off at any point (%, %), arrives at the point (--Uo, 
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—Vo) with a slope just equal and opposite to the slope at {u^y, I'o), or arrives 
at the original point (uq, Vo) with the original slope in the same direction, 
corresponds to terminal motion. Of course, a process of trial and error has 
to be adopted in order to obtain such a path. 

If, in addition to gravity, other forces passing through the centre of mass 
act on the lamina, the process is essentially the same. In the case of a 
propeller thrust, this additional force makes a constant angle with the 
direction of the lamina, so that in each interval the representative point in 
the (u, 'y)-plane travels an additional distance proportional to the thrust 
parallel to this constant direction. 

The method also applies to the square or nearly square lamina of small 
dimensions with moderate rotation. For a wider lamina, or a large square 
lamina, the method needs modification. 

It seems that some work on the wide lamina has been done by Gerlach, 
Ahlborn, and Koppen, see ' Encyk. d. math. Wiss.,' IV, 17, pp. 172-3. 

14. In figs. 7, 8 the path and direction of a narrow lamina of considerable 
aspect ratio moving under gravity have been plotted for different initial 
conditions. In each case the arrowheads represent the direction of the lamina, 
and the points at the centres of the arrowheads show the centre of mass of 
the lamina. Gravity is downwards, and the rotation is clockwise. In each 
case the motion is supposed to commence at the point where the letter 
designating the case is placed. The early part of the motion in G is shown 
enlarged in G'. 

The cases we have plotted are sufficient to indicate the nature of the 
motion under practically any initial conditions and for any length of time. 
The initial part of the motion depends upon the initial conditions. After some 
time, whose lengtli depends upon the velocity initially, the path approxi- 
mates to an undulating curve about a line drawn downwards, to the left if 
the rotation is clockwise, and to the rigiit if the rotation is anti-clockwise. 
This is the terminal motion for the given value of A/co. 

It will be noticed that for initial conditions intermediate between those of 
A and B of fig. 7 a terminal motion would be obtained at once. The subject 
of terminal motion is discussed elsewhere by the author for simplified law of 
resistance. 

15. All the paths in figs. 7, 8 have been obtained for the same values of the 
resistance coefficient A, and the angular velocity oo. But we can readily 
show that the nature of the motion is the same for all values of A and w, 
and that in order to obtain the paths for all possible values of A and cd, we 
have only to use different values of P2P4 on the standard scale of ^g. 4, 
bearing in mind in each case what the actual scale of this plane is in the 
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special case under consideration. When this is done, it is found that the 
actual ('2t, v) curves for motion under gravity retain the same general shapes 
as those from Which the paths in figs. 7, 8 have been plotted. The general 
motion under gravity is therefore always of the same nature, being charac- 
terised by an ultimate oscillation about a line inclined to the downw^ard 
vertical on the same side as the rotation. 

The general similarity between the fundamental curves of figs. 4 and 5 
shows that this general property is also true of the nearly square plane 
lamina of small dimensions. 

A few simple experiments with a light lamina, as, for instance, a thin 
postcard, wdll corroborate the results obtained by the graphical method. 
The characteristic oscillation about a line on the same side of the downward 

vertical as the rotation is observed in all 
cases of the two-dimensional motion of such 
a card given an initial spin. 

In order to obtain more exact verifica- 
tion of these conclusions, photographs were 
taken at very short intervals of a card 
moving in the manner discussed. In this 

., , . , , work the author w^as joined by Dr. A. M. 

: ■ 'i^v^r^/^^i^' ^'^ Tyndall, D.Sc, Head of the Physics Depart- 
ment of the University of Bristol. Several 
of the paths plotted in figs. 7, 8 were 
obtained. One of these photographs is 
reproduced in fig. 9. 

The case of a lamina^ with negligible 
moment of inertia, corresponding to Lan- 
chester's phugoid theory, and of a lamina 
oscillating about a vertical fall, correspond- 
ing to the case of the parachute, have been 
investigated, and will be published else- 
where. 
It is with much pleasure that I acknowledge my great indebtedness to 
Prof. Bryan, to whom the initiation of this work is due. Prof. Bryan not 
only suggested to me the problems discussed in this and other papers, but he 
also gave me the benefit of his advice and criticism throughout the whole of 
the investigations. 
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